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We explore the possibility of observing elliptic flow in low multiplicity events in central pp collisions 
at LHC energy, -^5=14 TeV. It is assumed that the initial interactions produces a number of hot 
spots. Hydrodynamical evolution of two or more hot spots can generate sufficiently large elliptic 
flow to be accessible experimentally in 4-th order cumulant analysis. 
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Recent experiments at RHIC produces convinc- 
ing evidences that a collective matter is created in 
Au+Au/Cu+Cu collisions [J-Q. The evidences come 
mainly from observing finite, large elliptic flow in non- 
central collisions. Elliptic flow is the 2nd harmonic in 
the Fourier expansion of the momentum distribution of 
the identified particles, 
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4>r being the azimuthal angle of the reaction plane. El- 
liptic flow (iaj) measure the azimuthal correlation of pro- 
duced particles with respect to the reaction plane. Finite 
elliptic flow is now regarded as a definitive signature of 
collective effect 0, E| . It is also best understood in a col- 
lective model like hydrodynamics Q. In a non-central 
collision, the reaction zone is spatially asymmetric. Dif- 
ferential pressure gradient convert the spatial asymme- 
try in to momentum asymmetry. In other words, in a 
hydrodynamic model, spatial asymmetry of the interac- 
tion region produces collective effects. Now protons also 
have finite extension (though of smaller size) as do a nu- 
cleus. It is then possible that in finite impact param- 
eter pp collisions at the Large Hadron Collider (LHC), 
at cm. energy y / s=14 TeV, asymmetric reaction zone 
will produce collective behavior which could be mani- 
fested as elliptic flow. However, even if flow is produced, 
whether or not it will be accessible experimentally will 
depend on both the flow strength and the multiplicity 
in the phase space window where the flow is measured. 
This is because non-flow effects like di-jet production, 
also show azimuthal correlation not related to the reac- 
tion plane. They need to be disentangled for faithful 
reconstruction of the reaction plane. Several standard 
methods have been devised to discriminate non- 

flow effects. Event plane method [f| Q determine the 
reaction plane, but require large multiplicity for unam- 
biguous determination. Cumulant method [9| does not 
require measurement of the reaction plane. Cumulants 
of multiparticle azimuthal correlation are related to flow 
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harmonics. The cumulants can be constructed in increas- 
ing order according to the number of particles that are 
azimuthally correlated. The method relies on the differ- 
ent multiplicity scaling property of the azimuthal corre- 
lation related to flow and non-flow effects. In the cumu- 
lant method, for charged particle multiplicity n mu i t , v-i 
can be reliably extracted using two particle correlator, if 
t>2{2} > 1/ \jn mu it- Higher order correlators will increase 

the sensitivity, e.g. t>2{4} > l/ 71 ™*;*- Maximum sensi- 
tivity «2 > could be achieved using still higher 
order correlators (cumulants of order greater than 4) . In 
the Lee- Yang zero method [ijj elliptic flow is obtained 
from the zeros in a complex plane of a generating func- 
tion of azimuthal correlation. It is also less biased by the 
non-flow correction, V2 > ^-/n mu u- It may be mentioned 
here that the sensitivity arguments are based on order 
of magnitude and are mostly valid in large multiplicity 
limit. 

Simulations of pp collisions at LHC energy by event 
generators like PYTHIA, indicate that while charged 
particle multiplicity in central rapidity region peaks at 
n m uit < 10, the distribution has a pronounced high mul- 
tiplicity tail. There are appreciable number of events 
with multiplicity n mu i t > 50, comparable to multiplicity 
in peripheral Au+Au collisions at RHIC. If high multi- 
plicity (n mu it > 50) pp collisions generate elliptic flow 
v 2 > 0.15, flow would be experimentally accessible in the 
2nd or higher order cumulant analysis. However, bulk of 
the pp collisions produces multiplicity n mu i t ~ 10. For 
n mu it ~ 10, one readily gets v 2 {2} > 0.3, v 2 {4} > 0.17, 
v 2 {> 4} > 0.1 as the scales of elliptic flow needed to be 
produced to be experimentally accessible with 2nd, 4th 
and higher order cumulant method. Note that the min- 
imum «2, V2 ~ 0.1 that can be accessed experimentally 
(with higher order cumulants) is rather large. Peripheral 
Au+Au collisions at RHIC energy produces much less 
flow, v 2 ~ 0.06 [Hi. 

In the present paper, in a hydrodynamic model, we ex- 
plore the possibility of observing rather large elliptic flow 
(v2 > 0.1) in low multiplicity, n m ult ~ 10, pp collisions 
at LHC. Applicability of hydrodynamics in a small sys- 
tem like pp is uncertain. Hydrodynamics require 'local' 
thermal equilibration, which can be achieved only if the 
mean free path of the constituents is small compared to 
the size of the system, A << R. In pp collisions, size 
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of the system is not large, A ~ R ~ 1 fm. However, it 
can be argued [l2| that the essential assumption of ideal 
hydrodynamics is that the stress tensor is isotropic in 
the local rest frame, Ty = pSij, with some equation of 
state relating pressure p to energy density. The condition 
Tij = pdij is a statement of isotropisation of the medium. 
If the medium is isotropised within a time scale Ti, hy- 
drodynamic may be applicable beyond Tj. In a QGP, 
non-abelian version of Weibel instabilities can grow very 
fast, isotropising the medium 12]. The maximum growth 
rate is, 7 ~ g\/n/p hard , n being the density and Phard is 
the characteristic momentum scale dominating the exci- 
tations in non-equilibrium QGP [12]. In the saturation 
scenario, Phard ~ Qs, the saturation scale [l3| . If one 
further assume, n ~ Ql/g 2 , then 7 ~ Q s , so the typi- 
cal isotropization time scale is 1/Q S . In a central (b=0) 
collision, this scale in proton is larger than in a nucleus, 



Qs,A > Qs,p- UM- Taking Q s ~ lGeV, one obtains the 
isotropization time scale as ~ 0.2 fm, justifying applli- 
cability of hydrodynamics in central pp collisions after 
^ « 0.2 fm. Indeed, similarities between pp and Au+Au 
collisions has been observed even at RHIC energy (l5j . 
When phase space restriction due to conservation laws is 
taken into account transverse momentum distribution in 
pp and Au+Au collisions at RHIC energy show similar 
behavior (l5j . However, similarity in px spectra alone 
does not prove that collective model like hydrodynamic 
is applicable in pp collisions. Observation of finite elliptic 
flow could be a definitive signature of collective behavior 
in pp collisions. Recently, several authors have 
applied ideal hydrod yna mics to study elliptic flow in pp 
collisions at LHC. In [l7l[l8| . it was conjectured that high 
multiplicity events will show large elliptic flow, v-i ~ 10- 
20% if hot spot like structures are produced in the initial 
collisions. In [l||, smooth initial energy density configu- 
ration was evolved hydrodynamically. For energy density 
appropriate for pp collisions at LHC, in peripheral colli- 
sion small elliptic flow ~ 2-3% , was predicted. With 
smooth energy density, hydrodynamical evolution donot 
generate elliptic flow in central collisions. 



Presently, we assume that in a 'central' pp collision 
at LHC energy, a number N s of 'hot spots' are created. 
Recently, in [13, [l8| , scenarios with hot spots formation 
are considered. We do not dwell on the mechanism of 
hot spots formation. In the constituent quark model, 
only a few partons (quarks) interact, forming flux tubes 
between them. The hot spots can be formed from the 
rapid decay of the flux tubes [HI- However, relevance of 
constituent quarks in LHC energy collisions is question- 
able. The collisions dynamics will be governed mainly 
by gluons. Large (gluon) density fluctuations can lead 
to formation of hot spots pjj ■ Hot spots are assumed to 
have Gaussian density distribution with cr=0.3 fm. For 
N s number of hot spots, the energy density of the system 
can be obtained as, 
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FIG. 1: (color online) Transverse profile of the energy density 
in a random event, with (a) one, (b) two, (c) three and (d) 
four hot spots, are shown. Peak energy density is assumed to 
be £o=100 GeV/fm 3 . 
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The centre of the hot spots (r^) can be anywhere in the 
reaction volume. Hydrodynamic evolution of the fluid 
will depend on the positions (r.;). In the following, we as- 
sume that (rj)'s are randomly distributed within a sphere 
of radius R=1.2 fm. In Fig. [1] transverse energy density 
profile of a random event, with one, two, three and four 
hot spots are shown in four panels. The peak energy den- 
sity is assumed to be Eq— 100 GeV/fm 3 . As it will be 
discussed below, for e =100 GeV/fm 3 , boost-invariant 
hydrodynamical evolution of 2-3 hot spots, produces 
charged particle multiplicity in the range n mu it = 8 — 10, 
close to the number obtained in simulations with event 
generators. 

From FigO] it is obvious that depending on the number 
of hot spots, initial spatial asymmetry or the participant 
eccentricity will vary. For a finite number of hot spots, 
participant eccentricity can be obtained as (l9j . 
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where at = 



; (* 2 )-(*) 2 , al = (x 2 )-(x) 2 , al = (y 2 )-(yr, 
a X y = (xy) — (x}(y) and (...) denote density weighted 
averaging. In the limit of homogeneous density distribu- 
tion, e coincides with the usual definition of participant 

eccentricity 6], e/ = Tpfefey- I n table jll we have noted 
the participant eccentricity as a function of the number 
of hot spots. They are averaged over 200 random ini- 
tial configurations. Error are statistical only. When only 
a single hot spot is formed, e is consistent with zero. 
Participant eccentricity is finite if 2 or more hot spots 
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FIG. 2: Probability distribution P(v 2 
with two, three and four hot spots are are shown in panels (a- 
c). P( v 2n A J l '^ lt ) in mixed events I and II are shown in panels 
(d) and (e). In panel (f), -P(«2W^* it ) in viscous (r//s=0.08] 
fluid evolution, with three hot spots, is shown. 



are formed. It also appear that for more than one hot 
spots, participant eccentricity is approximately constant, 
e ~ 0.5, which is comparable to that in very peripheral 
Au+Au collisions. 

We assume that 'baryon free' hot spots are 'locally' 
thermalised in the time scale r,=0.2 fm and evolve hy- 
drodynamically. The space time evolution of the fluid is 
then obtained by solving the energy-momentum conser- 
vation equation, 



(4) 



where T^ u — {e-\-p)u ,J 'u 1 ' —pg^ v , is the energy-momentum 
tensor, e, p and u being the energy density, pressure and 
fluid velocity respectively. Assuming boost-invariance, 
Eqs0]is solved in (r = \Jt 2 — z 2 , x, y, rj s = i In ) co- 
ordinates, with a code " 'AZHYDRO-KOLKATA" 2 ', de- 
veloped at the Cyclotron Centre, Kolkata. Initial fluid 
velocity was assumed to be zero. Details of the code 
can be found in [20l - |22j . Solution of Eq|4] requires an 



equation of state (EoS). We assume a lattice based EoS 
with connnement-deconfinement cross over transition at 
T co =196 MeV. Details of the EoS can be found in (HHH 
Recent lattice simulations [24| have been parameterised 
to obtain EoS for the deconfmed phase. EoS of the con- 
fined phase is that of a non-interacting hadronic reso- 
nance gas comprising resonances with mass m < 2.5 



TABLE I: Event averaged (for sample size of 200) participant 
eccentricity (e) are shown as a function of number of hot spots 
(7V S ). Also shown are the event averaged charged particle 
multiplicity ((n mll (t)), mean transverse momentum ((pt)) and 
elliptic flow (v2) as a function of number of hot spots (N s ). 
The bracketed quantities are the averages when sample size 
is halved. Last two rows corresponds to mixed events EI and 
mixed events EII (see text). 
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1 





4.97 ±0.02 


0.722 ±0.001 


0.003 ± 0.001 






4.97 ±0.02) 


(0.722 ±0.001) 


(0.003 ± 0.001) 
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0.532 ± 0.052 


7.75 ±1.17 


0.634 ± 0.054 


0.147 ±0.071 






(7.88 ± 1.11) 


(0.632 ± 0.054) 


0.152 ±0.068) 


3 


0.536 ±0.051 


9.68 ± 2.24 


0.599 ± 0.037 


0.160 ±0.053 






(9.87 ±2.12) 


(0.601 ± 0.040) 


(0.158 ±0.056) 
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0.457 ±0.048 


11.05 ±2.58 


0.582 ± 0.029 


0.161 ±0.050 






(11.39 ±2.67) 


(0.581 ± 0.026) 


(0.160 ±0.049) 


EI 




8.36 ±2.91 


0.634 ± 0.065 


0.118 ±0.019 


EII 




8.45 ± 2.36 


0.627 ±0.057 


0.138 ±0.022 



GeV. Hydrodynamic evolution of initial QGP fluid, with 
the lattice based EoS correctly reproduces large vol- 
ume of experimental data in Au±Au collisions at RHIC 

HQ. 

For a given hot spots configuration, hydrodynamic 
equations are solved to obtain freeze-out surface at a fixed 
temperature Tp=130 MeV. Using the Cooper-Frey pre- 
scription, invariant distribution of 7r~ is obtained. In the 
present paper, we have neglected resonance contribution. 
We multiply the ir~ multiplicity by a factor of 1.5 to ap- 
proximately account for the resonance contribution and 
noting that only ~80% of charged particles are pions, 
it is further multiplied by a factor of 2 x 1.2 to obtain 
the charged particle multiplicity. In tableU event aver- 
aged charged particle multiplicity (n mu it), mean trans- 
verse momentum (pt), and transverse momentum inte- 
grated elliptic flow (V2) arc noted as a function of the 
number of hot spots. Sample size is N event =200. The er- 
rors shown are statistical. In tableU we have also noted 
the averages when the event size is halved (the brack- 
eted values). It is worth noting that even if the sample 
size is halved, averages and statistical error remain ap- 
proximately unchanged. This indicates that the sample 
size is statistically large and simulation results are sta- 
ble. The general trend of the (n mu i t ), (pr), and (1*2) are 
well established. Multiplicity increases with number of 
hot spots N s , e.g. as N s increases from 1 to 4, (n mu it) 
increases by a factor of ~ 2. The result is understood. 
We have fixed the peak energy density £o=100 GeV/ fm 3 . 
As the number of hot spots increases, initial total energy 
of the system also increases, and so does the multiplicity. 
For N a =l, the average multiplicity (n mu u) ~ 4. Aver- 
age multiplicity is nearly a factor of 2 less than the peak 
multiplicity (~7-8) expected in pp collisions. However, 
multiplicity will increase if £q is increased and expected 
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charged particles multiplicity could be reproduced even 
with a single hot spot. In pp collisions, one expects 2-3 
hot spots are formed. For 2-3 hot spots with peak en- 
ergy density 100 GeV/ fm 3 , average number of charged 
particle multiplicity (n mu i t ) = 8 — 10 is close to the ex- 
pected multiplicity. While (n m uit) increases with number 
of hot spots, mean transverse momentum of charged par- 
ticles decreases. However, for two or more hot spots, the 
decrease is less than 10%. It appears that the slope of 
the transverse distribution of charged particles will re- 
main approximately constant if two or more hot spots 
are formed. Simulation results for elliptic flow are most 
interesting. Determination of elliptic flow requires the az- 
imuthal orientation of the reaction plane. For a number 
of hot spots, arbitrary located, the azimuthal orientation 
of the reaction plane is non-trivial. We calculate Vi as- 
suming that the reaction plane is oriented such that the 
vi is maximised. As expected, (vz) is consistent with 
zero if a single hot spot is produced in the collisions. 
For a single hot spot, participant eccentricity is approx- 
imately zero and hydrodynamic evolution does not gen- 
erate any flow. Indeed, events with a single hot spot 
essentially corresponds to smooth hydrodynamics, and 
with respect to elliptic flow development, are equivalent 
to events without any hot spot. If single hot spot forma- 
tion or smoothed hydrodynamics dominates, elliptic flow 
will not be observed in central pp collisions. If in pp colli- 
sions more than one hot spot is formed, then substantial 
flow is generated. It also appear that, for two or more hot 
spots, (V2) is approximately constant, (i^) ~ 0.15 ±0.05. 
With two or more hot spots in the initial state, central pp 
collisions at LHC could generate more elliptic flow than 
that generated in peripheral Au+Au collisions at RHIC. 

We have also considered mixed events, (i) EI, when 
probability of formation of one, two, three and four hot 
spots are equal and (ii) EII, when probability of forma- 
tion of one, two, three and four hot spots are 10%, 50%, 
30% and 10% respectively. In a realistic situation mixed 
events EII are more likely than EI. In table|H in last two 
rows, (nmuit) , (Pt)j and (v^) for the mixed events EI and 
EII are noted. In mixed events EI and EII, elliptic flow 
is still large, (v 2 ) ~ 10-15%. 

Simulation results indicate that if more than one hot 
spot is formed in initial pp collisions, large elliptic flow 
may result in low multiplicity events. However, generat- 
ing large elliptic flow does not ensure that they will be 
accessible in experiment. As discussed earlier, whether 
or not the flow is accessible experimentally depend on 
both the flow strength (i^) and multiplicity (n m uit)- For 
n mu it ~10, event plane method @, Q of determination of 
the reaction plane will have very large uncertainty. Cu- 
mulant method [9] do not measure the reaction plane. To 
find the order of cumulant required for faithful measure- 
ment of V2, we computed the measure (i) V2-y/n mu it and 

(ii) V2-T?Jy^i t . Simulated events donot satisfy the con- 
dition V2-y / n mu it > 1, necessary for two particle cor- 
relator to measure t>2- In Fig|5] probability distribu- 
tion P{v2-T\l^ lt ) of ^2-"muit m ^ ne simulated events are 



shown. In panels (a-c) distribution obtained with two, 

three and four hot spots are shown. P(v2-i\lf dt ) in mixed 
events I and II arc shown in panel (d) and (e). In all the 
cases, we find that more than 20% events crosses the 
threshold, f2 > 1 f° r 4th order cumulant analysis. 

In the present simulations, we have neglected the ef- 
fect of viscosity. Recently, from a systematic analysis of 
STAR data on <f> mesons, QGP viscosity was estimated, 
v/s = 0.07 ± 0.03 ± 0.14 [22], the first error is statistical 
and the second one is systematic. The central value of 
the estimate is close to the ADS/CFT lower bound on 
viscosity, n/s > 1/4tt (25|. Other parameters remaining 
unchanged, effect of viscosity is to reduce elliptic flow and 
increase particle multiplicity. For example, in Au+Au 
collisions, compared to ideal fluid evolution, in viscous 
fluid (rj/s = l/47r) evolution, elliptic flow is reduced by 
~10%, and particle multiplicity is increased by ~ 20% 
22]. However, if hot spots are formed in pp collisions, 
the gradients will be steeper and it is possible that vis- 
cous effects will be more than that in nuclear collisions. 
Indeed, explicit simulations also indicate that viscous ef- 
fects are more in pp collisions. For example, in a typical 
pp event, with three hot spots, with peak energy density, 
£o=100 GeV j fm 2 , compared to ideal fluid evolution, in 
minimally viscous (rj/s— 0.08) fluid evolution, elliptic flow 
is reduced by ^30%, and particle multiplicity is increased 
by ^80%. If one adjust the initial energy density to ob- 
tain multiplicity (n mu i t ~10), as in ideal fluid evolution, 
then elliptic flow is reduced by ~15%. In FigJ^f, we have 
shown the probability distribution, P{v2.n^t lt ), in mini- 
mally viscous (j]/s= 0.08) fluid evolution when three shot 
spots are formed in initial pp collisions. The peak en- 
ergy density is adjusted to £o=50 GeV/ fm 3 , such that 
average particle multiplicity is n mu it — 11.39 ± 0.83. 
Approximately 15% of the events satisfy the criterion, 
3 / 4 ^ 1 

In the present simulations, we have considered only 
hot spots with size er=0.3 fm. The size of the hot spots 
is important for the development of elliptic flow, and it is 
important to understand the effect on simulations if hot 
spot sizes are changed. In a hydrodynamic model, elliptic 
flow depends on the initial spatial eccentricity. Initial 
spatial eccentricity increases if hot spot size is reduced 
[17j . Elliptic flow in the model will increase if hot spot 
size is reduced. Explicit simulation indicate that elliptic 
flow will increase by ~10% if hot spot size is reduced 
from 0.3 to 0.2 fm. The criterion V2n Z J l ^ llt > 1 will be 
better observed then. On the contrary, if hot spot size is 
large, V2 decreases and the criterion V2T\j l ^ ilt > 1 will be 
less well observed. 

Before we summarise our results, we would like to com- 
ment on 'typicality' of hot spot formation. We have con- 
sidered central pp collisions. Central collisions do not 
dominate the cross section. If events with at least two 
hot spot are only a small sub set of already rare central 
events, then the present result, large elliptic flow in cen- 
tral pp collisions due to hot spot formation, may not be 
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relevant experimentally. We argue here that two hot spot 
formation probability is not small, rather it is large. Hot 
spot formation is possibly connected with multiple inter- 
actions of partons. For example, double partonic scat- 
tering (DPS) may lead to formation of two hot spots. In 
hadronic collisions, evidences of double parton scatter- 
ing have been found in multi jet events [111 Double 
parton scatering events are not rare. For example, CDF 
collaboration [27| estimated the double parton scattering 
events as ~ 52% in p + p — > 7 + 3jets + X reactions at 
v / s=1.8 TeV. Theoretical investigations also predict large 
number of double parton events at LHC energy [2^, [2£j ■ 
In minimum bias collisions, more that 50% double parton 
scattering are predicted. In events with hard process, the 
fraction is even larger. 

To summarise, due to finite extension of protons it is 
possible that at LHC energy (y/s—lA TeV), p p co llisions 
will produce collective matter. Recently in jl7l [l8j . it 



was conjectured that if pp collisions led to formation of 
hot spot like structures, in high multiplicity n mu it > 50 
events, collective effects will be manifested as 'experimen- 
tally measurable' elliptic flow. Bulk of the pp collisions 
are low multiplicity (n mu u ~ 10) events. We have ex- 
plored the possibility of observing elliptic flow in low mul- 
tiplicity events. In a hydrodynamic model, we have simu- 
lated pp collisions. It is assumed that initial collisions led 
to formation of a number of hot spots, distributed ran- 
domly in the interaction volume. Hydrodynamical evo- 
lution of two or more hot spots with peak energy density 
£o=100 GeV/ fm 3 , produces charged particle multiplic- 
ity, n mu i t rj 8 — 10 as expected in bulk of the events in 
pp collisions. Formation of two or more hot spots also 
generate substantial elliptic flow, V2 ~ 0.15 ±0.05. Large 
elliptic flow in low multiplicity pp collisions will be ac- 
cessible experimentally in 4th order cumulant analysis. 
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